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Abstract 

We analyze the sequence of polynomials defined by the differential- 
difference equation Pn+i{x) = Pn{x) + x{n + asymptotically 
as n ^ oo. The polynomials Pn{x) arise in the computation of higher 
derivatives of the inverse error function inverf(a:). We use singularity 
analysis and discrete versions of the WKB and ray methods and give 
numerical results showing the accuracy of our formulas. 
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1 Introduction 



The error function erf(x) is defined by pLj 

X 

erf(x) = —= I exp (— t^) dt 



(1) 



and its inverse inverf (x), which we will denote by satisfies 3 [erf(a;)] = 
erf [3f(a;)] = x. The function appears in several problems of heat con- 
duction |12j. In [10] we considered the function 



N{x) 



27r J -oo 

and its inverse S{x), satisfying 

S [N{x)] = N [S{x)] = X. 
It is clear from ([1]) and ([2]) that 



(2) 



erf (^1+1 



N{x) = 

S{x) = V23 {2x - 1) . 



and therefore 

In [To] we showed that 

S'{x) = V^TT exp 

and 



2 ^ ' 



(3) 

(4) 
(5) 



^(") = P„_i(S)(5')" n>l, 
where Pn{x) is a polynomial of degree n satisfying the recurrence 

Po(x) = l, P„+i(a;) = P^(x) + x(n + l)P„(x), n > I. (6) 

The same approach was employed by Carlitz in |2j . From ([6]) , it follows easily 
that for a fixed value of n 



Pry(x) n\ x"" , X ^ OO. 



(7) 
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From ([3]) and we conclude that 



= 2^P„_i (v^j) {3') 



n 



n> I. 



Since 



3{0) 



0, a'(o) 



a/tt 



2 



we have 




) 



'P„_i(0). 



(8) 



It follows from ([8]) that estimating 3^"'\0) for large values of n is equivalent to 
finding an asymptotic approximation of the polynomials Pn{x) when x = 0. 

The objective of this work is to study Pn{x) asymptotically as n ^ oo 
for various ranges of x. We shall obtain different asymptotic expansions for 



— > oo and (i) < x < oo, (ii) x = O {n ^) and (iii) x = O ^A/ln(n) j . The 
paper is organized as follows: in Section [2] we approach the problem using a 
singularity analysis of the generating function [T3] of the polynomials Pn{x)- 
In Section [3] we apply the WKB method to the differential-difference equation 
(El). In [15j, we used this approach in the asymptotic analysis of computer 
science problems and in [6] to study the Krawtchouk polynomials. Finally, 
in Section H] we analyze ([6]) again using the ray method [13| and obtain an 
asymptotic approximation valid in various regions of the (x, n) domain. In 
[1], 05 Hi we employed the same technique to analyze asymptotically other 
families of polynomials and in [S], [H] to study some queueing problems. 

2 Singularity analysis 

In [TDj we obtained the exponential generating function 





which implies that 
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where the integration contour is a small loop around the origin in the complex 
plane. Using (jl]), we have 



PJx) = e-^'l^^- 



27ri / 

|2|<r 



1 dz 



-a;2/2. 



\z\<r 

n\ / n + 1 



dS\N{x) + zN'{x)\ 



27ri / 

|2|<r 



■S\N{x) + zN\x)\dz 



and therefore 



(n + 1)! 
27ri 



5[iV(x)+ziV'(x)]- 



(9) 



|z|<r 



Since «S'(a;) has singularities at a; = and x = 1, we consider the functions 

N{x) 



1 - N{x) 



N' x) 



N'(x) 



We have 



and 



Zi(-x) 



v^exp - Zi 



[x) 



(10) 



= x-^ + 0{x-^) , 



X ^ 00, 



Changing variables to 



277 exp ^ + + O (a;~^) 



m; = [z- Zi{x)]N'{x) 



X ^ 00. 



in (IHD, we obtain 



P. fx) 



1 (n + 1)! 



A^'Cx) 27ri 



c 



w 
N'ix' 



+ ZAx) 



-(n+2) 



4 



or, 



J{w + l) 



where C is a small loop about w = w*{x) in the complex plane, with 



w (x) 



X). 



To expand f|TT]) for n — >■ oo with a fixed a; G (0, oo) , we employ singularity 
analysis. The function J{w) has singularities at w = ±1. By ([T]), we have 



J 

2 / _i2 , ^ _ -72 

W7 = —= e dt = 1- e ^ 

'■K J 




+ o {j-^) 



J ^ oo, 



so that 



J^{w) ~ a/— ln(l — w), w — i> 1 
and by symmetry we have 

J{w) ^ -^/^h^(TTw), w^-l+. 

The integrand in (ITT]) thus has singularities at w = and w = — 2, but 
for x > 0, the former is closer to w*{x). We expand (fTTl) around = by 
setting w = 5/n and using 



^i(x) 



77 ^2/2^ 

n 



-(n+2) 



[Zi{x)] 



-(n+2) 



exp 



^l(x) 



• (12) 



Then, we deform the contour C in (ITTI) to a new contour Ci that encircles 
the branch point at if = (see Figure [I]) . This leads to 

P„(x) ~ (n + l)!v^e^'/2 1 (13) 

00 r 2 " 

/" (T+-T-)exp 



where 



T±(5, n) = Viivr - In (S) + In (n). 
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Figure 1: A sketch of the contours C and Ci. 

Here T^{6,n) corresponds to the approximation of J^{w + 1) for w ^ 0, 
above or below the right branch cut in Figure [H 
For n large we have 



T+{6,n) -T-{6,n) 



TTl 



n] 



and then evaluating the elementary integral in (|T3l) leads to Pn{x) ~ '^i{x, n) 
as ^ oo with 



-(n+l) 



n+l 



(14) 
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Figure 2: A plot of In [P4o(x)/40!] (solid line) and In [^i(a;, 40)/40!] (ooo). 



and 




(15) 

In Figure E] we plot In [P4o(x)/40!] and In [^i(x, 40)/40!]. We see that the 
approximation is very good for x = 0(1) but it becomes less precise as 
X oo. This is because our previous analysis assumes that n ^ oo with 
< x < oo. If either a; or x ^ oo, we must modify it, which we will do 
next. 

When X 0, or more generally when x = O (n^^) , the singularities at 
w = and w = —2 are nearly equidistant from w*{x). On the scale x = y/n, 
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y = O (1) , we have 



Z,{x) = zJ^) = J--^ + 0{n-'), n^oo (16) 



n+l 

n\ I I2\ (12 



and (fT4|) simplifies to 



'■'^''^'^v?m[r.) ^H'v^'- '''' 

But, to this we must add the contribution from w = —2, which corresponds 
to replacing y by —y and multiplying by (—1)" the right hand side of (fT7|) . 
We note from ffTOl) that 



y|e^"/2«; + Zi(-x) = ^e-'/' (w + 2) - Z,{x), 



so that the integrand in (|TT|) is antisymmetric with respect to the map 
(x, w) — s> (—X, — 2 — w) . Thus, for n — s> oo and x = O {n^^) we get ~ 
\E'2(a;, n) with 



^2(1/, n) 



n+l 

n! / /2" 



^2 In (ra) IV TT 



exp(?/y^) + (-irexp(-i/y^ 



(18) 

As ?/ — s> 00, the alternating term becomes negligible and (ITSl) matches to (fT^ . 
as X 0. In Figure E] we plot the ratio In [P40 {^) /40!] / In [^2(1/, 40)/40!] 
and verify the accuracy of ( fTSl) . 

Letting x — > cxd in f|T3|) using f|T4l) yields 



Pnix) 



n] 



which differs from ([7]). This suggests that another scale must be analyzed, 
where x and n are both large. Thus, we consider the case of x — > 00, with 
X = O ^Vln 71^ . Now the singularity at w = in (ITTl) becomes close to w*(x), 
since 

*( \ ^ 
w [x) ' 



Xx rn 
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Setting z = ^/x with ^ = 0(1), we obtain 

S[N{x) + zN\x)] = V2J [l + xl^ 

= V2J 



i + o(.-)) 




~ V2 J Y + ln(i) - i 1.1 0) - 1.1 (1 - f ). 

Thus, we have 



Pn{x) 



27ri 



Ci 



'x^ , X 

h In 

2 Vl-e 



Here the contour Ci is a small loop about = 0. Now we again employ sin- 
gularity analysis, with the branch point at ^ = 1 determining the asymptotic 
behavior for ri oo. A deformation similar to that in Figure [H leads to 



Pnix) 



n\x 



1 



V2 + ln(nx)-iln(^) 



(20) 



For X >> A/ln(n) this collapses to (JTj) . 

By examining (HM and (l20l) . we can obtain the following approximation 



PJx) ~ ^■:^(x,n) 



+ 2 In(nx) — In (^^^ 



n n+l 



(21) 



which is more uniform in x, since it holds both for x = 0(1) and x = 

O ^Vln for n large and for x — > oo with n fixed. However, we must still 

use (1181) if n is large and x is small. In Figure H] we plot In [P4o(x)/40!] and 
In [\E'3(x, 40)/40!] and confirm that fl2T|) is a better approximation than ([1] 
for large values of x. 
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80^ 






60^ 












20^ 









/ 2 4 


6 8 10 

X 


Figure 


4: A plot of In [P4o(x)/40!; 


(solid line) and In [^3(0;, 40)/40!] (000). 



3 WKB analysis 

We shall now rederive the results in the previous section by using only the 
recurrence and ([7]). We apply the WKB method to IQ, seeking solutions 
of the form Pn{x) = n\Pn{x), with 

Pn{x) ^ exp[{n + l) A{x)]B{x,n), n ^ 00. (22) 

Thus, we are assuming an exponential dependence on n and an additional 
weaker (e.g., algebraic) dependence that arises from the function B{x,n). 
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Using ([22]) in dSD leads to 



d 



B{x,n) + —B{x,n) 
on 



O 



1 d 

= b + A'ix)] B(x, n) H TrBix, n) 

n + lox 

Expecting that || = o (5) and = o (|f ) , we set 

e^(-) =x + A{x) 

and 



on n + lox 



^ B(x,n) ~ —-^B(x,n) 
n ox 



(23) 
(24) 



To solve (!23l) we let 



to find that 



X 



A{x) = — — + a{x) 



a'(x) = e-"'/2e«W. 



Solving this separable ODE leads to 



X 



A{x) = -—-\n [C {x) + k] 



(25) 



where is a constant of integration. To fix k, we assume that expansion 
(122|) . as X ^ oo, will asymptotically match to ([7j), when this is expanded for 
n — s> oo. In view of ( [221) this implies that 

Pn{x) ~ = exp [?T,ln(x)] , X ^ oo, 

so that 

74(x) ~ ln(x), X oo. 
In view of f l25|) this is possible only if A; = and then from f|T5|) we have 



A(a;) 



In 



(x) ~ ln(a;), x ^ oo. 



(26) 



We next analyze ( j24l) . Using ( !26l) to compute 6"^*^^^ we obtain 

-r^X7^5(a;,n) = -— 5(x,n). 
C (x) on n ox 



12 



Solving this first order PDE by the method of characteristics, we obtain 



B(x, n) 



n 



where 6 (■) is at this point an arbitrary function. However, since n is large and 
X = 0(1), we need only the behavior of 6 (■) for large values of its argument. 
We again argue that by matching to ([7]) we have 

exp [{n + 1) A{x)] B{x, n) ~ x", a; — > oo, 

and using fl26l) we get 

1 



Bix, n) ~ e 



-A(x) 



X ^ oo 



X 



and thus 



so that 



b (^nxe^ 
biz] 



X oo 



X 



Combining our results, we have found that 



2; — > oo. 



PJx) 



V2^1n(n) - In [C (x) 



ax) 



n+l 



n 00. 



(27) 



This applies for x = 0(1) and n ^ 00, where we can regain the results of 
the singularity analysis (IT^ by simply using 



A/ln(n) — In [( (x)] ~ \/ln{n), n 



00. 



Formula (1271) is also valid for n = 0(1) and x — > 00, where it reduces to ([7]). 
However, (!22|) breaks down for x ^ 0. 

We thus consider the scale x = y /n, y = O (1) and set 



P„(x) = n\Pn{nx) = n\Pn{y), 



with which ([6]) becomes 



Pn+i [y + - 

n 



n + l n 



(28) 



(29) 
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For fixed we seek an asymptotic solution of fl29p in tlie form 

^n(l/)~e""g(i/,n), n oo, (30) 



where q{y,n) will have a weaker (e.g., algebraic or logarithmic) dependence 
on n. From fl2^ we obtain, using fl3Ul) . 



yd 

Qiv, n) + T-q{y, n) + -^q{y, n) + (n^^) 
on nay 



(31) 



1 - - + O (n-2) 



-K-q{y,n) + -q{y,n). 
oy n 



If g (y, n) has an algebraic dependence on n, then ^q{y, n) should be roughly 
O {n~^) relative to q{y,n) , ^q{y,n) roughly O {n~'^) and so on. Thus, we 
expand q {y, n) as 



q {y, n) = qo{y, n) + -qi{y, n) + O {n 
n 



(32) 



where qo{y,n),qi{y,n) have a very weak (e.g., logarithmic) dependence on n 
and balance terms in (|3T|) of order 0(1) and O (n~^) to obtain 



and 



d 

e''qo{y,n) = —qo{y,n) 



qi{y,n) + -^qiiy,n) + y-^qoiy,n) 
^?i(y, n) - -^qoiy, n) - yqo{y, n). 



Solving (l33l) yields 



qoiy,n) = exp (e"y) q(n) 



(33) 



(34) 



(35) 



where q(n) must be determined. We could solve (IMI) . using (!35|) . but its 
solution would involve another arbitrary function of n. Thus, considering 
higher order terms will not help in determining (]{n). Instead, we employ 
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asymptotic matching to (I27|) . Expanding fl27p for x ^ and comparing the 
result to (1251) as ?/ — > oo, with (15U]) . (15^ and (1531) . we conclude that 

« = iiiifzV qH= J-^ - (36) 



2 V TT 



But then our approximation for y = 0(1) is not consistent with P„(0) = 
= P„(0) for odd n. We return to (l29l) and observe that the equation also 
admits an asymptotic solution of the form 

P„(l/)~(-l)"e^"g(l/,n), n ^ oo 

where, analogously to (!33l) . we find that 

so that another asymptotic solution to (l29l) is 

P,(y)~(-l)'^e^"exp(e^y)q(n). (37) 

We argue that any linear combination of (l30l) and (|37|) is also a solution and 
that the combination which vanishes at j/ = for odd n has /3 = a and 
q(n) = q(?T.), as in (|36|) . We have thus obtained, for y = 0(1), 



p.. n! /2 



a/vt ln(n) 




exp \y\ -] + (-l)"exp 




n ^ oo. 



This agrees with ( fTSl) . obtained by singularity analysis in section 2. 

To summarize, we have shown how to infer the asymptotics of Pn{x) using 
only the recursion ([6]) and the large x behavior (JTj). Our analysis does need 
to make some assumptions about the forms of various expansions and the 
asymptotic matching between different scales. 



4 The discrete ray method 

We shall now find a uniform asymptotic approximation for Pn{x) using a 
discrete form of the ray method fllj. This approximation will apply for x 
and/or n large. We seek an approximate solution for ([6]) of the form 

Pnix) = exp[f{x,n) + g{x,n)], (38) 
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where g = o{f) as n — > oo. Since P„(x) = x", n = 0, 1 we see that we must 
have 

f{x,n) ~ n\n{x) and g{x,n) (39) 
as n ^ 0. Using (138|) in ([6]), we have 



fdf Id^f dg\ fdf dg\ , , 
as n ^ oo, where we have used 

df 1 

/(x, n + 1) = f{x, n) + —{x, n) + 2 ^(^' ^) • 

From ( l40l) we obtain, to leading order, the eikonal equation 

| + („+l).-exp(f)=0, (41) 

and the transport equation 

Id^f dg dg f df\ ^ 

+ TT" - Tf - exp --^ = 0. (42) 



2 c?n^ dn dx \ dn^ 

To solve fHTl) . we use the method of characteristics, which we briefly 
review. Given the first order partial differential equation 

F(x,n,/,p,g) = 0, with P=^, ^=1^' 

ox on 

we search for a solution /(x, n) by solving the system of "characteristic equa- 
tions" 

dx dF dn dF 

dt dp ' dt (9g ' 

dp _ dF OF dq _ OF dF 

dt dx ^ df dt dn ^ df 

df _ dF dF 

dt ^ dp ^ dq^ 

with initial conditions 

F [x(0, s), n(0, s), /(O, s),p{0, s), g(0, s)] = 0, (43) 
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and 



^/(O, s) = p(0, s)^x(0, s) + g(0, s)^n{0, s), 



(44) 



where we now consider {x, n, f,p, q} to all be functions of the variables t and 
s. 

For the eikonal equation fHT]) . we have 



F (x, n, f,p,q) = p — e'^ + {n + 1) X 
and therefore the characteristic equations are 



dx ^ dn 
dt ' dt 



dp / , dQ 
— = -n + l), — 



(it 



dt 



and 



(it 



p — qe* 



Solving (146|) subject to the initial conditions 

a;(0,s) = s, n(0,s) = 0, p(0,s) = A(s), g(0,s)=5(s) 
we obtain 



x = t + s, n = -W^exp ( y + 5 



erf 



t + s 



V2 



erf 



^expf^ + 5)(t + s) 



1 



erf 



t + s 



From ( !39|) we have 

A{s) = and 5(s) = ln(s) 
which is consistent with (l43l) . Therefore, 



+ exp ( --f -st + B] - t-e^ + A, 



-h^ -st + B 



x = t + s, ra = -^/^sexp(y 



erf 



^^-erf^^ 



V2 



(45) 

(46) 
(47) 



(48) 



(49) 



(50) 



p 



|sexp f Y ) (t + s) 



e,.f I ^i^) - erf (-^= 



+sexp ( --t^ - st ) - (t + s) , g 



-it2-st + ln(s). 
2 



(51) 



17 



Figure 5: A plot of the rays x(t, s),n(t, s) for s G [—2. .2] 



In Figure [5] we sketch the rays x{t, s),n{t, s) for s G [—2. .2] . 
Using (ICT) in (H7I) we have 



df pn f \ , ^ 



+ s 



l + -t^ + st- In (s) 



exp ( --t^ -St] -{t + s). 



(52) 



Using (jlH]) in (ilj), we get 



/(0,s) = /o 



(53) 
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and solving (152]) subject to (1551) . we obtain 



f{t,s) 



exp , 
2 ^ V 2 



erf 



X s 



or, using ([SD, 



1 + -t^ + - In (s) 



1 



- ( 2^' + 1 + ^ 



/ = [In {s)-l]n + - {s^ - x^) {n + l) + U 



To solve the transport equation (H2|) . we need to compute | 
as functions of t and s. Use of the chain rule gives 



and hence, 



" 8x 


ax ' 




'9t_ 


av 




'1 


0" 








ax 


an 










ds 


as 







1 


.at 


as. 




.ax 


an. 





-at 


dt_- 


1 


an 


ax- 




a,n 








.ax 


as 
an. 


~ J{t,s) 


an 
. at 


ax 
at . 



where the Jacobian J{t, s) is defined by 

dx dn dx dn dn dn 
^' ^ dt'd^ ~ 'd^'dt ^'d^~~dt' 
Using we can show after some algebra that 

J = i s -\ — jn + s. 



Using q = ^ in fj42l) . we have 



2 dn dn dx 



or 



d_ 

dn I 2^' 



di _ dg_^g 
dx dn 



and using (H6l) . we obtain 
d (\ 



dn\2 



dg dx dg dn dg 
dx dt dn dt dt 
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Since — = we have 



dn 



2 
1 

2J 



1 d f dn 
~ ^2d^ \ dt 
d'^n dx d'^n dx 



1 / d'^n dt d'^n ds 
+ 



+ 



dt'^ ds dtds dt 
1 d f dn 



2 V dt'^ dn dtds dn 
d'^n 

+ 





1 / d^n 


•)= 


2J V dt^ 


dn\ 


1 dJ 


'dtj 


~ ~2J~dt' 



dtds 



where we have used (!56l) and (1571) . Thus, 

dg _ 1 dJ 
di ~ ~2J~dt 



and therefore 



g(t,s) = --HJ)+C{s) 



for some function C{s). Since from fl39|) we have (7(0, s) = 0, while 
J(0, s) = s, we conclude that C(s) = | ln(s) and hence 



Using fl58|) we can write (|59|) as 



2 



J(t,s)J 



2 



(ra + 1) + n 



Replacing / and g in fl38|) by fl55|) and fl6Q|) . we obtain P„(a;) ~ $ (a;, 
n — cxD, with 



Kexp 



1 



(s^ - x^) (n + 1) - 



where k = e^" is still to be determined. 
Eliminating t from ([5]) we get 



n + 



s exp 



erf 



71; 



+ 1 + 



- 
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which defines s{x, n) imphcitly. For every n > there exist only two solutions 
Sm{x,n) < and Sp{x,n) > of (1^ (see Figure). Since erf(x) is an odd 
function, it follows that 



-Sp{-x,n) 



(63) 



Although we have P„(x) ~ ^ [x,n; Sm{x,n)] for x ^ —1 and P„(x) ~ 
^ [x,n; Sp{x,n)] for x ^ 1, the two approximations are comparable when 
X is small and therefore we must add their contributions. 

We shall now find the constant k, in (1611) by using ([7]). We rewrite (I62p as 



' exp (s^) 



n 



-\ -2 



exp 



n2 



(64) 



and for a fixed value of n, consider tlie limit x ^ oc. It follows from (E 
tliat Sp {x, n) ~ X and therefore we consider an expansion of the form 



Sp{x, n) X + So + six + S2X ^ + S3X~ 



X 



00. 



(65) 



Using fl65|) in fl64|) . we obtain 
■So = 

S3 = 1 — In (ra + 1) — 
and therefore 



S2 = 0, si = In (n + 1) , 
In^ (n + l) 



s^exp (s^) ~ (n + l)'x2e"' 



X 



n 



00. 



Solving (|66|) we have 



5'p(x, n) ~ Y W [(n + 1) x^e^^], x — > cxo, 
where W (z) denotes the Lambert W function, defined by [3] 

W {z) exp [W (z)] =2, e C 
and having the asymptotic behavior 

In In(z)"' ^ 



(66) 



(67) 



W iz) = ln(z) - In In (z) + + Q 

m{z) 



00. (68) 
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Using (1^ and (1^ in (IM]) . we obtain 



$ (x, n; s) ~ K (n + l)""^^ g '^3;'^^ x — > 00. 



From Stirling's formula 



V27m + O (^n 5^ 



n"e n ^ cx), 



(69) 



and 

+ 1)"+ 

we conclude that 
and thus 

$ {x, n; s) = exp 



e\/n + O in 2 



(s^-a;2-2) {n+1] 



(ra + 1) + n 



Using fl63|) . we have Pn{x) ~ ^4(0;, n) as n ^ 00, with 

\l/4(x, n) = $ [x, n; Sp{x, n)] + $ [x, n; —Sp{—x, n)] , n ^ 00. 



(70) 



(71) 



In Figure [6] we compare ln[P4(x)/4!] and In [\&4(x, 4)/4!] for < x < 10 and 
in Figure [7] for — 1 < x < 1. We note that the asymptotic approximation flTTl) 
is more uniform than f|T^ . ( |T8|) and ( l20l) but it is less explicit since Sp{x,n) 
must be obtained numerically. 

Next, we compare the results of this section with those in the previous 
two sections. We first consider x > 0, with x = 0(1) and n — >• cx). From (162|) . 
we have 



Sp{x,n) 



\ 



W 



jn + l) 
e(x) 



(72) 



where ({x) was defined in (ITSl) . Using (1721) and (l68l) in (!70l) . we get 

,n n+l 



Pn(x; 



riTT 



In ( 72 ) 



C(:^) 



which agrees with (HM after taking flB^ into account. 
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Figure 6: A plot of In [P4(x)/4!] (solid line) and In [^4(0;, 4)/4!] (000). 

Now we consider the limit n —>■ 00, with x = y/n and y = 0(1). From 
fl62|) . we have 



S,{y/n,n) ~ ^/W (—) + (l + \^y] • (73) 

V ^ ; V ^ vr y (^)n 

Using ([75]), dHSD and dHSD in dZOD, we find that 



n 

2n I I2\ I 12 



*(»/«, n;S„) ~ I I exp I -J-y 



ln(n) I V TT y I V tt 



23 




Figure 7: A plot of In [P4(x)/4!] (solid line) and In [^4(0;, 4)/4!] (000). 
and therefore 

agreeing with f|T8l) . 

Finally, we consider the limit n ^ 00 with x = u^y\n{n), u = 0(1), 
u > 0. From fl62l) . we have 

Sp (u^/\n{n), ~ [(n+ l)^^^"' ln(n)] . (74) 
Using dZH) and dMl) in dZO]), we have 

Pn yu^yin {n)j ~ n''e~''V27m-^=== y^/ln {n)j , n ^ 00, 
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which agrees with fl20l) . 
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